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Abstract 

We propose a model based on the primitive system of the Navier-Stokes equations in 
a bidimensional framework as the I - plane approximation, which allows us to explain 
the variety of tracks of tropical cyclones (typhoons). Our idea is to construct special 
analytical solutions with a linear velocity profile for the Navier-Stokes systems. The 
evidence of the structure of linear velocity near the center of vortex can be proven 
by the observational data. We study solutions with the linear-velocity property for 
both barotropic and baroclinic cases and show that they follow the same equations 
in describing the trajectories of the typhoon eye at the equilibrium state (that relates 
to the conservative phase of the typhoon dynamics). Moreover, at the equilibrium 
state, the trajectories can be viewed as a superposition of two circular motions: one 
has period 2n/l, the other one has period 2n/bo, where I is the Coriolis parameter 
and bo is the height- averaged vorticity at the center of cyclone. 

Also, we compare our theoretical trajectories based on initial conditions from the 
flow with tracks obtained from the observational database. It is worth to mention 
that under certain conditions our results are still compatible with observational 
data although we did not truly consider the influence of steering effect. Finally, 
we propose the parameter-adopting method so that one could correct the weather 
prediction in real time. Examples of our analysis and the use of parameter-adopting 
method for the historic trajectories are provided. 
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1 Introduction 



Tropical cyclones, commonly known as hurricanes in the North Atlantic Ocean 
and typhoons in the western North Pacific Ocean, are one of the most dev- 
astating weather phenomena in the world. The intense winds associated with 
tropical cyclone often generate ocean waves and heavy rain, which usually re- 
sult in severe disasters. To reduce the harm caused by the tropical cyclones, 
the most important issue is the prediction of its moving trajectory. 

For a long time, the phenomenon of formation of a huge intense atmospheric 
vortex attracted attentions of scientists working in various fields (e.g. see [1] 
and references in [2]). A large progress has been made in the last decade in 
understanding the physical process and the quality of operational prediction 
of tropical cyclones ([3], [I], and references therein). In particular, statistical 
and numerical weather prediction models have been improved dramatically 
and used to provide operational forecasters the high-quality prediction for 
national weather bureaus. 

The system of equations served for the numerical prediction is very compli- 
cated. It needs to take many factors into account, for examples, phase tran- 
sitions, vertical advection, and boundary layer effects, etc. Generally, it is 
difficult to use such system to comprehend its underlying physical processes. 
To have a conceptual understanding about these processes, some assumptions 
must be made to simplify the equations. The solutions then may appear in a 
simpler form so that it becomes much easier to gain their physical interpreta- 
tions. However, neglecting some terms in the equations may risk losing some 
important properties of the system, e.g. the symmetries. Therefore, it is very 
important to find particular solutions for the untruncated system of equations 
and to adjust them to track some particular phenomena. In the present paper, 
we seek special solutions for the primitive air-motion system. As we should 
show later, these solutions can be considered as a model of the typhoon when 
it is near its center. 

Our model is bidimensional and we does not consider the effect of heating that 
results from a coupling between the latent heat released in the clouds and the 
vertical wind shear. It is known that a great variety of the vortex phenomenon 
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can be explained in the frame of two-dimensional models if we consider the 
vertical dimension as a small scale comparing with the horizontal ones (e.g. 
[5]). In meteorology, this concept is often adopted in studying the processes 
for middle and large scales p]. In this case, the atmosphere is regarded to 
be "thin". For example, the meteorologists often consider the atmosphere as 
a one-layer fluid in the barotropic framework without specifying the vertical 
coordinate. This bidimensional theory, despite such simplicity, can use to ex- 
plain many observed features of the tropical cyclone motions [3] . However, the 
essence of bidimensionality in our model differs from above. We derived it from 
a primitive complete three-dimensional system of equations by averaging the 
system over the height, so the vertical processes of the system are presented 
in a hidden form. 

The paper is organized as follows. In section 2, we study the derivation of 
a bidimensional model of the atmosphere dynamics based on the primitive 
Navier- Stokes equations for compressible viscous heat-conductive gas in the 
physical 3D space. In section 3, we investigate polynomial solutions for both 
barotropic and baroclinic system and find that the equations which govern the 
trajectories at the equilibrium point are the same in both cases. In section 4, we 
test possible trajectories of the typhoon center with various input parameters 
and demonstrate that our the artificial trajectories can imitate the behavior of 
the typhoon eye tracks. Moreover, numerical simulation shows that the motion 
of the vortex center provided by our explicit algebraic formulas is reasonably 
close to the motion of a developed vortex when it is near the equilibrium. 
In section 5, we examine two different types of historical typhoons: Man Yi 
(parabolic trajectory) and Parma (loop), with our theoretical results. Here, 
we propose a method of the parameter-fitting that allows us to predict the 
trajectory of typhoon eye only based on three successive points in the past. In 
section 6, we further investigate the situation when the typhoon encounters 
the shore. In this case, the typhoon may disappear due to the destruction of 
the stable vortex. Finally, discussions about our topics and future work are 
provided in section 7. 



2 Bidimensional models of the atmosphere dynamics 

The motion of compressible rotating, viscous, heat-conductive, Newtonian 
polytropic gas in R x M 3 , is governed by the compressible Navier-Stokes equa- 
tions [7] 

d t p + div(pu) = 0, (2.1) 
d t (pu) + Div(pw ®m) + V?) + 2wx« = DivT, (2.2) 

dt Q/°M 2 + P e ) + div [i\p\ u \ 2 + pe + p)u^j = div(Tw) + V(«, V0), (2.3) 
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where p, u = (ui, U2, U3), p, e, 9 denote the density, velocity, pressure, internal 
energy and absolute temperature, respectively. Here u = (cji, u; 2 , cj 3 ) is the 
angular velocity of the Earth rotation, k = (/c 1( k 2 , k 3 ) 7 K{ > 0, i = 1, 2, 3, is 
the vector of heat conduction. The stress tensor T is given by the Newton 
law 

T = Tij = p (diUj + djUi) + X divu 5ij, (2.4) 

where the constants p and A are the coefficient of viscosity, and the second 
coefficient of viscosity. Here Div and div stand for the divergency of tensor 
and vector, respectively, <S> denotes the tensor product of vectors. 

The state equations are 

p = Rp9, e = c9, p = Aexp(^\ p 1 . (2.5) 

Here A > is a constant, R is the universal gas constant, s = log e— (7— 1) log p 
is the specific entropy, c = -J^, 7 > 1 is the specific heat ratio. 

The state equations (2.5) imply 

p = (7 - l)pe, (2.6) 

which allows us to consider (2.1) - (2.3) as a system of equations for the 
unknowns p, u, p. Indeed, from (2.1) - (2.3) and (2.6), we arrive 

d t p + (u, Vp) + ipdivu = (7 - 1) E TijdjUi + £ A-. (2.7) 

i,j=i R P 

For simplicity, we denote the system (2.1), (2.2), (2.7) as (NS). 

Let x = (xi,X2) be a point on the Earth surface, <^o be the latitude of 
some fixed point x . Following from [8] (early this approach was used in [9] 
for barotropic atmosphere), we derive a spatially two-dimensional system for 
(NS). Let us introduce cf> and / to represent for taking the average of and / 

over the height respectively, namely, <p := / (f)dz, f :— — / pf dz. where 

Jo p Jo 

(p and / are arbitrary functions, and denote g{t, x) = p, P(t, x) = p, U(t, x) = 

u, 0(t,x) = 9. Moreover, the usual adiabatic exponent, 7 is related to the 

27-I 

"two-dimensional" adiabatic exponent 7 as follows: 7 = < 7. 

7 

Now, we include the impenetrability conditions in our model. These conditions 
ensure that the derivatives of the velocity equal to zero on the Earth surface 
and a sufficiently rapid decay for all thermodynamic quantities as the vertical 
coordinate z approaches to infinity. In other words, the impenetrability condi- 
tions make sure the boundedness of the mass, energy, and momentum in the 
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air column. They also provide the necessary conditions for the convergence of 
integrals . 



Let us denote additionally I = 2u 3 sin (p Q (the Coriolis parameter) and L = 

-1 

1 

In this way, one can get the system as the I- plane approximation near x : 
g(d t V + (U ■ V)U + ILXJ + VP) = Ti(U), 

d t Q + V • (f?U) = 0, (2.8) 



d t P + (U, VP) + 7 PdivU = ( 7 - 1) T 2 (U) + «(4 A6 - 0, 

ti 



where 



Ti (U) = DivTjj, T„ = /i + 9,-Ui) + A divU 

T 2 (U) = ((2// + AXO^U!) 2 + (<9 2 U 2 ) 2 ) + M<9iU 2 + ^Ui) 2 ), 

fj, and A are the coefficients of viscosity interpreted as the coefficients of tur- 
bulent viscosity (may be not only constants). Their values are much greater 
than those for the molecular analogues. For simplicity, the coefficients of the 
heat conductivity are assumed to be constants and equal. £ is the heat flow 
from the ocean surface, it is also assumed to be a constant. 



3 Models of the typhoon dynamics with the linear velocity profile 



Suppose that the velocity vector near the origin has the form 

u(t,x) = a(t)r + b(t)r ± , (3.1) 

where 

r = (xi,x 2 ) T , r ± = (x 2 , -xif. 

This assumption is made according to the observational data [16] (see FigJTJ. 
Moreover, as it is shown in [T2], in the stable vortex, the velocity field neces- 
sarily satisfies the Cauchy-Riemann conditions, namely, 

dx\ dx 2 ' dx 2 dx\ 
This property holds for (13. ip and also verifies our choice for the velocity field. 

Note that if we do not take into account the turbulent viscosity (/i = A = 0) 
and the heat transfer (k = 0) in our model. We get the system similar to the 
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Figure 1. The tangential component of (relative) velocity in m ■ s 1 on different 
stage of the typhoon formation (August 30, August 31, September 2). The hurricane 
Anita, North Atlantic, 1977 [E] 

gas dynamics instead of (2.1)- (2.3). Solutions to such system were constructed 
and also studied in [T3].[T1].[T5]. 

We next consider two forms of the system (2.8), which have exact polynomial 
solutions. 



3. 1 System with respect to the velocity, pressure and entropy, without viscos- 
ity and heat conductivity 

7-1 

We assume /i = A = k = and introduce a new variable it = p n . We also 
consider the 2D entropy S connecting with q and P through the state equation 
(2.5), where we use 7 instead of 7. For the new unknown variables 7r, U, S, we 
arrive the system [2] 

d t V + (U • V)U + ILU + exp -Vvr = q^TAU), 

7 - 1 7 

d t n + (Vtt • U) + (7 - 1)tt divU = 0, 
|(c^+(U,V£))=T 2 (U). 

Following (TOUT], we change the coordinate system so that the origin of the 
new system is located at the center of the typhoon eye. Now U = u + V, where 
V(t) = (Vi(t), V^i)) is the velocity of the typhoon eye propagation. Thus, we 
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obtain a new system 

d t u + (u • V)u + V + lL(u + V) + — ^— exp -Vvr = e _1 T 2 (u), (3.2) 

7 - 1 7 

<9i7r + Vvr • u + (7 - l)vrdivu = 0, (3.3) 

|(«9 t S+(u,VS))=T 2 (u). (3.4) 

Note that in [TUffTT] . shallow water equations were considered. It is a particular 
case of (2.8) with 7 = 2, S = const., A = [M — K — 0. 

Given the vector V, the trajectory can be found by integrating the system 

xi(t) = Vi(t), x 2 (t) = V 2 (t). (3.5) 



We assume the velocity has a linear profile (13. ip and look for other components 
of the solution to the system fl3.2l) - fl3.4p in the following form 

7r(t, x) = A{t)x\ + B{t)x 1 x 2 + C(t)x 2 2 + M(t)xx + N(t)x 2 + K(t), (3.6) 

00 

^(t,x) = S , (t)+ 7 ln J2 S ili2 {t)x[ix%. (3.7) 

U,i2=l 

Here, Ti(u) = ^(u) = 0. Recall that K(t) is the value of "renormalized" 
pressure n at the center of vortex, therefore, in the physical sense, K(t) > 0. 
Substituting fl3.ip - fl3.7p into fl3.2p - fl3.4p and match the coefficients for similar 
term, we obtain S ili2 (t) = 0, i%,i 2 G N, S (t) = const, A(t) = C(t), and 
B(t) = 0. At the center of typhoon, there exists a domain of lower pressure. 
Therefore, it is natural to set A(t) > 0. 

Note that the motion near the typhoon center is "barotropic" in our solution, 
i.e., the pressure only depends on its density. However, this barotropicity is 
only for the bidimensional density and pressure. This assumption may not 
hold for three-dimensional cases. 

Let us introduce a constant c := ^-exp 

Thus, the functions a(t), b(t), A(t), M{t), N{t), K{t), Vi(t), V 2 (t) satisfy the fol- 
lowing system of ODEs: 

A + 2^aA = 0, (3.8) 

a + a 2 - b 2 + lb + 2c A = 0, (3.9) 

b + 2ab-la = 0, (3.10) 

^ + 2(7-1)0^ = 0, (3.11) 

M + (27 — l)aM — bN = 0, (3.12) 

N + (2 7 - l)aN + bM = 0, (3.13) 
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Vi - IV 2 + c M = 0, (3.14) 
V 2 + lV 1 + c N = 0. (3.15) 

From (13. 8p and (13.101) . we have 

b= l - + C l \A\ ll \ (3-16) 

where C\. a constant. Therefore, system (13.81) - (13.101) can be reduced to 
equations 

72 

A = -2jaA, d = -a 2 - - + C\A 2h - 2c A. (3.17) 



Further, if A(t) and a(t) are known, we can find other components of solutions. 
Namely, from (I3.12p and (13.131) . c we can get 

M{t) = (M 2 (0)+iV 2 (0)) 1/2 exp ( J a(r)dr\ sin^t+d J* A 1 ^{r)dr+C 2 ), 







N(t) = (Af 2 (0)+iV 2 (0)) 1/2 exp / a(r)drj cos(^+C 1 f A 1 ^ (r)dr+C 2 ). 

From (13. 8p and (13.111) . we obtain 
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K(t)=C 3 (\A(t)\)- 

Here C 2 , C3 are constants depending only on initial data. However, from (13. 14ft . 
(I3.15p . and (13. 5p . the trajectory does not depend on K{t). 

The phase curves of (I3.17P can be found explicitly. They satisfy the algebraic 
equation 

a> = C^-CtAU 1 ^- T)+ ^ l A, (3.18) 

with the constant C4 depending only on initial data. Our analysis shows that 
there exists a unique equilibrium point on the phase plane (A, a), A > for 
1 < 7 < 2. It is a center on the axis a = 0. We denote this stable equilibrium 
as (A , 0), where A is a positive root of the equation — l -^ + C 2 A 2 ^ — 2cqA = 0. 



3.2 System with respect to the velocity, density and temperature 



We take the average of the first state equations (2.5) over the height to get 
P = RgO and then exclude the pressure from system (2.8), we arrive 

g(d t XJ + (U ■ V)U + ILU) + R(qV0 + tfV^) = Ti(U), 
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d t g+(Vg-U) + gVU = 0, 

a t e + (ve ■ u) + ( 7 - 1) e divu = ^— -t 2 (u) + ^(a6-0- 

K Kg 
After changing the coordinates system, it takes the form 

g(d t u + (u ■ V)u + V + lT(u + V)) + R(gV6 + 6Vg) = T x (u), (3.19) 

d t g+ (Vg-u) + gdivu = 0, (3.20) 

d t e + (ve ■ u) + ( 7 - 1) e divu = ^^t 2 (u) + -^(A0 - £). (3.21) 

The position of the vortex center can be obtained again from ( 13.51) . 

Note that now the viscosity and heat conductivity are not neglected. 

We can obtain a closed system of ordinary differential equations if we seek the 
solution to (I3.19I) - (I3.21I) in the following form: 

6(t, x) = A x (t)x\ + B x (t)x x x 2 + C x {t)x 2 2 + M 1 (t)x 1 + N x {t)x 2 + K x {t), (3.22) 

g(t,x)=K 2 (t). (3.23) 

Substituting (J3HD, ([3T22]) . f !3T2^|) into ( 13~T9l -( 13T2T]) . we obtain A x (t) = C x {t), 
B x (t) = and 

A x + 2 7 aAi = 0, (3.24) 

a + a 2 - b 2 + lb + 4RA X = 0, (3.25) 

b + 2ab-la = 0, (3.26) 

K 2 + 2aK 2 = 0, (3.27) 
k x + 2 {l -l)aK x = 2(7-1)(2^ + A) q2+ K (e _ 4A0; 

Jt A 2 

Mi + (2 7 - l)aMi - 6JVi = 0, (3.29) 

N x + (2 7 - l)aN x + bM x = 0, (3.30) 

V x - IV 2 + 2RM X = 0, (3.31) 

V 2 + IV X + 2RN X = 0. (3.32) 

Equations (13.241) . (13.251) . (I3.26[) can be reduced to the system of equations as 
before: 

72 

A x = -2^aA x , a = -a 2 -- + C 2 x A 2 x h - ARA X , (3.33) 
b= l - + C x \A x \ l l\ (3.34) 

where C x is a constant. 

Assuming that we know a(t) and A x (t), we can find b(t), M x (t), N x (t), V x (t) 
and V 2 (t) as before. The only stable equilibrium on the phase plane (a, A x ) 
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for 7 G (1,2) is a center (0,A ), where A is a root of the equation — — + 
<7?A2 /7 - 4iL4 2 = 0. 

It is important that the viscosity, the heat conduction and heat inflow terms 
show only in the equation (I3.28p . therefore, their values only affect the function 
Ki(t). Nevertheless, the shape of trajectory defined by the equations (13.241) - 
QSZSID , (15351) - ( 13321) does not depend on K x (t). 



4 Possible trajectories 



Let us analyze trajectories of a stable typhoon for both "barotropic" and 
"baroclinic" cases. At the equilibrium point, a(t) = 0. Other coordinates of 
the stable point for the "barotropic" system (p£Hj) . (j3~l?j) . fl37TU]) are A(t) = A 

and b(t) = bo = — + Ci(A )~ . In the real meteorological situations, I >> b . 
Thus, from (13TT21 - IXTSj) and (j3T5]) . we obtain 

, v V 2 (0) c M(0) . 
Xl (t)= Xl (0) + ^-l + ^l (4.1) 

(VM c N(0) \ . fV 2 (0) c M(0) \ 

+ l^-/(^J sm/t -l^ + /(^7)J cos ^ 

c A^(0) . , c M(0) 
6 (o - /) b (b - Z) 

„(,)= Xl(0) -m + «m (4.2) 

c M(0) . c N(0) 

■ sm b t + — — — cos b t, 



h(b -l) h(b -l) 



for I 7^ bo. The analogous results we get in the "baroclinic" case at the stable 
equilibrium for the system ( 13. 24ft . f l3.25[> . and (I3.26p . namely, a(t) = 0, Ai(t) = 

Ao, b(t) = bo = — + Ci(Aq)~ . To obtain the equations for the trajectories, we 
should change c M(0) to 2M 2 (0) and c N{0) to 2RN 2 {0). 

Thus, the trajectory is a superposition of two circular motions: one of them 
has a period 2tv/1, the other one has a period 27r/6 - They can lead to the 
appearance of loops, sudden change of directions and other complicated tra- 
jectories. Several examples of artificial trajectories are presented in Figsj2]and 
® 
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Figure 2. Possible artificial trajectories at 22°iV, 6 = 10 s , 
c = 10~ 1 ,M(0) = N(0) = 10 -9 , V2(0) = 0(in corresponding units) [2] 



V^OJ --0.63m Is, 




Figure 3. Possible artificial trajectories at 22° N, bo = 10 5 s 
c = 10" 1 ,M(0) = N(0) = 10~ 9 (in corresponding units) [2] 



It is possible to compute the position (xi,x 2 ) of the trajectory based on the 
system (13.81 - 1X11)1) . (I3.5p by applying ODE integrators. However, the difference 
between the numerical trajectory and points on the curve (4.1), (4.2) for the 
physical values of parameters is very little. For example, for the initial data 
a(0) = 10- 5 s-\ 6(0) = 5-10- 5 s-\ A(0) = 10" 9 (N-m.- 2 )^- 1 ^ -m" 2 , N(0) = 
10~ 3 (N ■ m~ 2 )^~^^ ■ mr 1 , M(0) = 2iV(0), K(0) = -1 ■ m ■ s~\ V 2 (0) = 
1 ■ m ■ s -1 , Xi(0) = x 2 (0) = 0, the difference in positions after three days is 
about 20 km (see Figf4]). (Here, the radius of typhoon is taken as 300 km, 
the difference between the pressure at the center of typhoon and the ambient 
pressure is 10 4 Pa, the constant c is estimated as 10 _1 ). Therefore, basically, 
there is no reason to get the solution curves by applying numerical integration 
for (13.81 I5TT3]) and (13.51) instead of using explicit formulas (4.1) and (4.2). 
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Figure 4. Difference between trajectories corresponding to the exact equilibrium 
point (solid line) and in a vicinity of equilibrium (dashed line) within three days 

5 Examples of analysis of real trajectories 

Assume initial conditions are known, one can integrate the system (13. SA - 
AB. 15p . (13. 5j) . But these initial conditions include the height-averaged physi- 
cal quantities such as divergency a(0), vorticity 6(0), parameters A(0), M(0), 
iV(0), (relating to the pressure field), and the initial velocity V(0). It is not 
clear how to measure these quantities. The only widely available information 
is positions of the typhoon eye. Thus, we try to fit the parameters using the 
formulas (4.1) and (4.2). However, it needs to assume that the vortex is at its 
conservative phase. 

We use the following algorithm: at the first step, we choose three successive 
points of trajectory ((aq (t ), a: 2 (to)), (^i(^i), ^2(^1)), (^1(^2), ^2(^2)))- Formulas 
(4.1), (4.2) give a system of four linear algebraic equations with respect to 
Vi(0), V^O), coM(O), coiV(O). The parameter 6 is free. The value of the Cori- 
olis parameter I is calculated at the initial point (xi(to), £ 2 (£o))- The second 
step is to choose the parameter b . If we deal with historical trajectory, b can 
be fitted so that the artificial trajectory coincides with the real one as long as 
possible. As we should show, the coincidence sometime lasts for several days. 

In practice, we can compute bo, for example, from the condition 



where ^(0) = Sil^gM , i = 1,2. However, ([EI]) and Q give different 
values of 60 basically. 



V 1 (0) = V 1 (0), 



(5.1) 



V 2 (0) = v 2 (o), 



(5.2) 
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Figure 5. The Man Yi trajectory [18] 



Let us estimate all possible values of parameter bo at the equilibrium point. 
As it follows from (13. 9p . b is a root of the quadratic equation 

b 2 Q - lb - 2c A = 0, 

therefore, bo ~ — 2c °, . for cqAq « I 2 . If we further assume A = 10~ 9 , cq = 
10 _1 , I = 10~ 4 , bo ~ —2 ■ 10~ 6 . Thus, to be physically meaningful, we can 
assume | £> 1 < 10 -5 and bear in mind that this value is essentially negative. 

Thus, if the set of three points is not able provide an appropriate value of b , 
the prediction is considered as a failure. We need to shift these three points 
until we get a definite value of bo that can be used for the weather forecast. 

We consider three successive points as suitable candidates for a forecast if 
equations (15. ip and (15. 2p have roots &oi and 602 in the interval |&o| < 10~ 5 
respectively and | £> i — &02 1 < £, with e sufficiently small. In this case, we chose 
bo as the mean of b i and 6 2- 

In the next, we give two examples of analysis of the real tropical typhoon 
tracks. 



5.1 Man Yi 



Let us consider a recent typhoon of the West Pacific region, Man Yi (4 cat- 
egory, July 8-15, 2007). The historic trajectory and observational data are 
given on Figs. 5 and 6 [18]. Remark that the numeration cited in Fig. 6 in- 
cludes certain intermediate points (like 1A, 2 A, etc). Here, We did numerate 
the intermediate points as in Fig. 6, instead, we adopt successive numeration 
in our examples. 
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Figure 6. The Man Yi data [U] 



3 5E6 



3EB 



; EES 



2ES 



1.5EB 



1EB 




1.45E7 1.5E7 1.55E7 1.SE7 



Figure 7. The Man Yi trajectory (observational data) with analyzed regions 
5.1.1 Historical trajectory: fitting averaged vorticity 



By analyzing the trajectory of the typhoon Man Yi, we found that it is possible 
to fit bo so that there is a good coincidence with the artificial trajectory given 
by formulas (4.1) and (4.2) in three different regions (FigJTj). Here and below 
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graphs, the diamond shape points represent the observational path of tropical 
typhoon and solid lines represent the analytic results; the first three points of 
trajectories, where the points and line coincide, are used for the parameters 
estimation. The first region corresponds to the first 14 points on FigJTJ (the 
total duration is 48 hours, the "tropical storm" stage). FigJH] shows the pre- 
diction for point 1 to 9 (18 hours), Figj9]- for point 1 to 14 (48 hours). Here, 
we get bo = 2.236198023 • 10~ 6 The second region on FigJT] corresponds to 
points 12 - 19 (33 hours, "tropical storm - typhoon 1", Fig JTO]) . At this region, 
the vortex is developing, thus we can not consider it as a stable one. Generally, 
making predictions according formulas (4.1), (4.2) may not give good results 
for unstable regions. Nevertheless, these formulas with bo = —1.8 • 10~ 5 s _1 
still gives a satisfactory coincidence for the real and artificial trajectories. The 
third region that we analyze is from point 22 to 34, which lasts for 72 hours. 
In this case, the typhoon is developing from 3th to 4th stage and then decays. 
We get b = —I - 10~ 5 and obtain a very good coincidence for the 3 day period. 




1 .61 E71.615E71 .62E71 .625E71 .63E7 1 .635E71 .64E7 



Figure 8. The first 9 points (18 hours). 
Here and below dots corresponds to ob- 
servational data, solid line represents an- 
alytic results 
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Figure 10. 12-19 points (33 hours) 




1.54E7 1.56E7 1.58E7 1.6E7 1.62E7 1.64E7 

Figure 9. The first 14 points (48 hours) 
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Figure 11. 22-34 points (72 hours) 
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5.1.2 Attempt of a forecast in real time 



Now, suppose that that we do not know the track of the whole trajectory in 
advance and try to find 60 to compare its value found from (15. ip and (15. 2p . 
It is worth to note that one may propose other methods to find the averaged 
vorticity. This will be a subject of our further investigations. Let us mention 
that one can consider the neural network approaches for this direction [T7]. 

If we set e = 2 • 10~ 6 , we get only two appropriate three point sets: one of 
them begins from the 21th point and gives a very good prediction (FigJTT]). 
here 60 ~ 1 • 10~ 5 s~ 5 ). The other one predicts a loop beginning from the 
point 7 (see Fig JT2|) . Thus, one can assume that at the 9th point, the vortex 
had gone through an exterior force. It constrains the trajectory to change the 
direction of its path. 

If we weaken the requirement for accuracy and set e = 3 ■ 10 -6 , we get one 
more three point set, beginning from point 24. In the latter case, the artificial 
trajectory gives a true direction. But comparing with the exact track, the 
artificial trajectory outgoing from the 21th point (Fig JTB"]) is bad. However, 
the spread of positions is inside of the annual mean diameter of the 33 percent 
strike probability region for 36 hours [I]. 





1.42E7 1.44E7 1 .46E7 148E7 1.5E7 1 .52E7 



145E7 1.5E7 1.55E7 1.6E7 



Figure 12. e = 2 • 10 6 : potential loop 



Figure 13. e = 3 • 10 : two best predic 
tions 



5. 2 Parma 



Now we consider the example of typhoon with a looping trajectory, Parma, 
20-31 October, 2003, FigJH FigED 

Fig. [16] presents the comparison of real trajectory within 78 hours of the ty- 
phoon in its conservative phase. The parameter b = — 6.1116738183-10 _5 s _1 is 
fitted for all points of trajectory according the historical data. Fig. [T71 presents 
a forecast for a long period (144 h.) based on the data for 12 hours. We can 
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Figure 14. The Parma trajectory |18j 
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Figure 15. The Parma data [18] 
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see that the forecast gives the true qualitative behavior of the trajectory, the 
loop. This loop was deformed by a steering flow. 




1.7E7 1.75E7 1.8E7 1.85E7 1.9E7 1.65E7 17E7 1.75E7 1.8E7 1.85E7 

Figure 16. The Parma trajectory from Figure 17. Forecast for 144 hours based on 
10 to 25th point (78 hours) 10-12 points (12 hours)) 



6 The role of surface friction 

It is well known that the typhoons basically do not last for a long time over a 
dry land. The key point is the significant increasing of the dry friction when 
the typhoon goes to the land. Now, let us add the damping term — kXJ in the 
right hand side of first equation (2.8), where A; is a nonnegative function of 
coordinates. For simplicity, we assume that A; is a constant. Therefore instead 
of ((335 JED), we get 

a + a 2 - b 2 + lb + 2c A = -ka, (6.1) 

b + 2ab - la = -kb, (6.2) 
equations (I3.8p . (l3.1ip - ( 13 . 1 5[) do not change. 

System (13. 8p . (16. ip . (16. 2p is closed. However, it does not encounter any equi- 
librium for i / 0. Therefore, it is not possbile to find a stable domain of 
low pressure from this system. Fig. 18 and 19 demonstrate the break-up of the 
stable equilibrium for k > 0. Computer simulations are made for initial data 
A = lCr 9 , a(0) = 0, 6(0) = -2 • lCr 6 , I = lCr 4 , c = 0.1 (in respective units). 
Fig. 18 presents the phase portrait of the system (I3.8p . (I6.ip . and (I6.2p for 
k — 0. In this stable equilibrium exists on the phase plane. Fig. 19 

shows the collapse of vortex for k > : the vortex strehgth slightly intensifies, 
whereas the motion becomes significantly convergent eventually. Here, initial 
data are the same as in Fig. 18, k = 3 • 1CP 5 . 



Thus, if we wish to obtain a stable vortex along with the surface friction, we 
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Figure 18. The phase portrait of (|3.8j) 
452) >OI for k = 0. 



Figure 19. The same an in Fig. 18 for 
k = 3 • 10~ 5 , computations are made for 3 
days. 



have to consider some additional phenomenon such as vertical advection (see 
[19]) in this context). 

Note that the system (13. 8ft . ( 16.11) . ( 16. 2 p is very interesting from the mathemat- 
ical point of view. It has a complex equilibrium at the origin and there is a 
possibility of existence of a strange attractor for certain values of parameters. 



7 Discussion 

In this paper, we show that the behavior of the developed tropical cyclone 
trajectory is determinated by exterior and inner parameters; the trajectory is 
a superposition of two circular motions: one has period 27r/Z, the other one 
has period 27r/6 - 

Our previous conclusion that the trajectories for both barotropic and baro- 
clinic models are governed by the same system of coefficient-equations seems 
conflict with the fact that baroclinicity has a strong relation with the cyclone 
formation. About this point, we would like to clarify that the cyclone forma- 
tion in our paper only refers to the conservative (or near conservative) phase 
of the atmosphere vortex development. Indeed, the evidence about the veloc- 
ity in a linear form ( 13. ip can be confirmed by the observational data which is 
collected as the developed typhoon is close to its center. As it shows in (Fig.l), 
at the stage of the vortex formation, the velocity may not to be linear. There- 
fore, we claim that we find a "toy" solution for a very complicated system to 
describe the processes of formation and decay of cyclones for which an analyt- 
ical solution is barely to be obtained. It worth to remark again that the term 
"barotropicity" in our paper is only restricted to describe the bidimensional 
framework, i.e. the term "barotropicity" is used with a hidden meaning of 
"taking an average over the height". In the physical 3D space, the flow does 
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not have to be barotropic. 



Our arguments can be verified from other researchers, for examples, in [20] . 
it says "The curvature of storm track is determined by a set of 'controlling 
parameters' related to the storm characteristics and ambient atmospheric cir- 
culation. These controlling parameters include the speed of storm (V(0) in our 
notation), storm intensity (6(0)), size (A(0)), the Coriolis parameter (/), am- 
bient atmospheric pressure field (M(0),iV(0)) and surface friction (k)" . Also 
in [3], it says that "In a barotropic framework, a tropical cyclone is basically 
"steered" by the surrounding flow (V(0), M(0), N(0)), but its movement is 
modified by the Coriolis force and the horizontal vorticity gradient of the sur- 
rounding flow." 



Certainly, it would be naive to claim that our model provides a better weather 
forecast for the typhoon trajectory than those modern models using numerical 
simulations. Firstly, the formula is derived under the assumption that the vor- 
tex is at stable phase, but, this assumption hardy holds for the real weather 
forecast. Therefore, we only can deal with oscillations near the equilibrium 
point at most in our model. Secondly, the only variables that we feed in the 
equations of thermodynamic parameters are initial data from real cases. How- 
ever, as the vortex moves, it may experience a forcing of baric fields, which 
may skew its trajectory significantly. Thus, generally, we can only expect to 
get the reproduction of the trajectories qualitatively, i.e., to predict a turning 
of the track without indicating the exact position of the vortex. Finally, we 
do not take into consideration the curvilinear geometry of the Earth surface. 
Therefore, only in a vicinity of an initial point, the linearization does not re- 
sult in a big distortion. The inclusion of the - effect does not give an exact 
solution with respect to the stable vortex (for approximative solutions, see 

!)• 



However, as in [3], even high-quality weather forecast models may not be re- 
liable due to its uncertainty about the initial conditions or for some other 
unknown reasons. Thus, operational weather forecasters still need to judge 
whether or not the prediction results should be taken into consideration to 
build an official weather forecast in their daily work. Therefore, it is helpful 
to have a complimentary tool to aid other existing models for the typhoon 
trajectory prediction. Our model seems very simple comparing with models 
used for the numerical weather forecast. However, it is very useful since it 
can be used to explain the trajectory behavior, predict the direction of the 
trajectory and the loop formation as well. 
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Our models can be further refined in two aspects: doing analytical refinement 
of the model and searching a better way for parameters-fitting. We conjecture 
that it is possible to obtain some solutions with analytical format for the 
Navier-Stokes system in spherical coordinates. Then we could include the /3 
- drift in our model and release the requirement that the solutions need to 
be bound around the neighborhood of the center of typhoon. The natural 
question arising in the parameters-fitting method is: can we actually use the 
average vorticity b as a measurement of predictability in the model? As it 
is mentioned in the previous sections, even feeding with the real data for 
the average vorticity b , one can only get approximate predication for the 
trajectory behavior. The quality of the predication strongly depends on the 
ambient meteorological fields, which are responsible for the steering effect. In 
particular, the "regularity" of the ambient meteorological fields seems to be a 
very crucial factor for the trajectory tracing. However, it can not be obtained 
only through analyzing the equations of trajectories. Thus, to improve the 
predictability and accuracy for our model, we need to combine our method 
with the analysis of available meteorological data and include other theorems 
like statistical moments (see in this context |22j) and neural networks |23j . 
[T7] in our model. 
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